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A variational model is used to calculate the static stability limits and 
equilibrium properties of "half -bubbles," magnetic domains residing on 
one surface of a magnetic bubble material platelet. Stability is achieved 
through the presence of gradients in the domain wall energy density and/ or 
saturation magnetic moment. The model evidences two distinct types of 
instability behavior separated by a critical value of the wall energy density 
gradient. Unlike the standard cylindrical domain, the half-bubble has a 
minimum stable value of the ratio of domain diameter to height. 

The half -bubble is shown to possess a number of properties which make 
it potentially useful for device applications. It is self -biased by its closure 
wall and in some cases is stable in zero external bias. Bias margins are of 
the same order as those for the standard cylindrical domain. It stabilizes 
on only one platelet surface, and its properties are independent of both 
material thickness and of minor irregularities on the second surface. In 
addition, its structure may be advantageous in avoiding the undesirable 
properties of hard bubbles. 

I. INTRODUCTION 

The current interest in cylindrical magnetic domains and their device 
applications has stimulated the search for such domains in a wide 
variety of materials produced by a number of growth techniques. In 
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some samples, Bobeck 1 has observed domains of reversed magnetization 
which apparently do not penetrate through the entire thickness of the 
sample (Fig. 1). These domains, referred to as "half-bubbles," have 
four features which make them potentially attractive for device applica- 
tions. First, they contact only one surface of the platelet in which they 
reside so that their properties should not be as critically dependent on 
surface preparation as those of the usual cylindrical domains. Second, 
as shown in this paper, their properties are independent of platelet 
thickness. Third, this article also shows that their top closure wall 
produces a self-biasing effect which makes external biasing unnecessary 
in some cases. Finally, Bobeck et al. 2 have observed in two-layer films 
that domains having a closure wall do not exhibit the undesirable 
properties of hard bubbles, and it is a reasonable presumption that this 
property also pertains to half -bubbles. 

Bobeck 1 has suggested that half-bubbles may be stabilized by the 
presence of gradients normal to the material surface in one or more 
of the material characteristics (e.g., domain wall energy density <r w 
and/or saturation magnetic moment M„). In this paper we use a varia- 
tional model to analyze the static stability of half-bubbles in materials 
having such gradients. Our model for the half-bubble and its instability 
modes is introduced in Section II. In Section III, we examine the stabil- 
ity of the model when the gradient in M , is zero. While no stability 
occurs in this simple case, the calculation of Section III is a useful 
preliminary to the more general calculation of Section IV, where the 
gradient in M, is assumed to be non-zero. The results of the stability 
calculations and their interpretation are discussed in Section V. 




Fig. 1 — Half-bubble magnetic domain configuration. 
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II. HALF-BUBBLE MODEL 

In constructing our model we idealize the half-bubble in two ways. 
First, since the details of the shape of the half-bubble are not known, 
we approximate the actual shape of the domain by a right circular 
cylinder of radius r„ and height h (h assumed to be less than the platelet 
thickness), as shown in Fig. 2. This approximation gives the model 
domain a form identical to that of the standard cylindrical domain, the 
stability of which has been calculated by Thiele. 3 We shall exploit 
this identity in calculating the half-bubble energy and its derivatives. 
Second, we defer all consideration of the microscopic structure of the 
domain wall of the half-bubble and assume that the wall energy density 
<t w varies linearly with distance through the platelet. 

The magnetostatic energy of our half-bubble model possesses an 
important invariance property. As we show in Appendix A for the 
case in which the M, gradient is zero, the magnetostatic energy of a 
cylindrical domain of fixed size located on the bottom surface of a 
platelet is independent of the position of the top surface of the platelet. 
This invariance is equally valid if M ', is a function of z, or if the domain 
assumes a more general shape. Thus, if we conceptually move the upper 
platelet surface down to the top of the model half-bubble in Fig. 2, 
we see that the magnetostatic energy of the half-bubble is equal to 
that of a standard cylindrical domain of the same size. 

This invariance property has two important consequences. First, in 
the limit where the M , gradient is zero, it allows us to obtain the 
magnetostatic energy of our model half-bubble from the results of 
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Fig. 2 — Model for the half -bubble magnetic domain. 
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Thiele 3 and Thiele et al. 4 for the standard cylindrical domain. Second, 
it indicates that the platelet thickness is not a fundamental length of 
the half-bubble as it is for the standard cylindrical domain. 

A convenient alternative standard of length, easily defined in terms 
of material characteristics, is provided by the material length at the 
bottom surface of the platelet (z = 0) : 

Z = <r w (0)/4*M*(0)- (1) 

In defining the a w and M , gradient parameters y w and y M , we then 
scale the variations in wall energy density and saturation magnetic 
moment to this material length: 

«r.« = ^(0)(l + y -f ) , (2a) 

M.(z) = M.(0)(l+^)- (2b) 

Experimental observation 1 indicates that there are three modes of 
instability which are of importance for the half-bubble: collapse, 
run-out, and "run-through". The general character of the collapse and 
run-out modes is already familiar from Thiele's stability calculation 
for the standard cylindrical domain; the half -bubble collapse is not 
purely radial, however, but involves axial collapse as well. The term 
"run-through" is used to describe the transition from a half -bubble to 
a standard cylindrical domain which contacts both surfaces of the 
platelet. 

As indicated in Fig. 2, the size of our model half-bubble is specified 
by the radius r and the height h. In considering stability against run-out, 
we include an elliptic distortion r 2 and specify the radial boundary of 
the domain by 

n(<p) = r + r 2 cos 2<p. (3) 

The model will be in equilibrium if the first derivatives of the total 
energy vanish. In the absence of in-plane anisotropy, dE T /dr 2 = at 
r 2 = by symmetry (cf. Appendix B). The equilibrium conditions 
are then 

dE T n dE T _ .. ,-. 

ir. = °' «r = °- (4a - b) 

The stability of the resultant equilibrium is determined from the second 
derivatives, stability occurring when 
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d*E, 

'd 2 j<: 



> 0, (6) 



and 



dK 



(*). 



> 0. (7) 



If the inequality in eq. (5) is replaced by the corresponding equality, 
then the half-bubble is on the verge of either collapse or run-through. 
If {dE T /dr\) eti , = 0, then the model is about to undergo run-out. 

III. HALF-BUBBLE STABILITY WITH CONSTANT M, 

We begin the actual calculation of the stability limits for the half- 
bubble by considering the limiting case of constant M, . This case 
lacks practical interest since we show that it leads to no region of 
stability for half-bubbles. In this limit, however, the problem loses much 
of the analytic complexity which characterizes the general case in which 
the gradient parameter y M j^ 0. Our method of solution-treating the 
second derivatives of the energy as polynomials in h/l with coefficients 
depending only on 2r //i-thus becomes clearly visible here. 

As noted in the preceding section, our half-bubble model is markedly 
similar in form to the model used by Thiele 3 and Thiele et al. 4 to treat 
the standard cylindrical domain. For this reason, in the present calcula- 
tion we adopt the notation of these papers wherever practical. 

The total energy of the half-bubble, measured relative to that of 
the uniformly-magnetized platelet, is 

E T = E w + AE„ + AE M (8) 

where E w is the wall energy, (including the energy of the top closure 
wall), AE H is the energy of interaction with the externally applied bias 
field, and AE M is the magnetostatic energy. Since the half -bubble is 
independent of the magnetic material above it, we can describe our 
model in terms of the unit step function u as 

M = $,[1 - 2u{r b {v) - r)]u(z)u(h - z) (9) 

where r b (<p) is as defined in eq. (3). 

The wall energy, consisting of integrals of a w (z) over the side and top 
surfaces of the half-bubble, assumes the form 
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E m = jT *.{z) dz fj \rl(<p) + (0j* dtp + wJk) fj ^ dp. (10) 

From eq. (10) it is clear that different wall energy densities for side and 
top walls could easily be included in our model. While such a distinction 
has intuitive appeal, experimental evidence on this point is lacking. 
We thus omit this added degree of freedom in the present calculation. 
The energy of interaction with the external bias field H = i,H, 
relative to that for the uniformly magnetized plate, is 

AE H = -2 f M-HdV = M.Hh f r 2 b (<p) dtp. (11) 

•'h.-b. »0 

If the elliptical distortion r 2 is set equal to zero, E w and AE H become 

E w = 27rr o ^cr„(0)(l + **) + in*r„(0)(l + ^) (12a) 

= 2rr haM2) + *rfr.(A), (12b) 

AE B = 2rr 2 hM.H. (13) 

It should be noted that the interaction energy AE H is proportional to 
the domain volume m*h. This characteristic is shared by the gradient 
term in the energy of the top wall [eq. (12)], leading us to view this 
term as arising from an effective bias field. The non-gradient part of 
the top wall energy lacks any /i-dependence, being proportional to irr 2 . 
It may be thought of as coming from an /i-dependent bias field analogous 
to that introduced by Liu 5 in connection with sputtered thin films which 
are exchange-coupled to a magnetic bubble material. 

The relative magnetostatic energy is obtained from the general 
expression 

B„-Udv(dV <*-*X*:; ih (14) 

4 J v J v > \T T J 

in the form 

AE M = M\ [ rdr [ d<p [ r' dr' 

Jq Jo Jo 

<k'[kk f - utp- 1 - (p 2 + h 2 y h ] (is) 



/ 



where 



k(r, r b (<p)) = 1 - 2u(r b (v) - r), (16) 
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p 2 = r 2 + r' 2 - 2rr' cos (<p - <p'). (17) 

This integral is given in closed form by Thiele et al. 3 for the case r 2 = 0: 
£E M = -2irh 3 -2irM 2 .I(x), (18) 

(19) 

where x = 2r /h and K, E are complete elliptic integrals. 

It is now a straightforward matter to differentiate E T with respect 
to r , h, and r 2 . Using the definitions for the normalized bias field S 
and Thiele's 3 wall-force function F 

= H/±tM. , (20) 

»» - ^ - ¥ K 1 + ?)X A) - • < 21 > 

we may conveniently write the equilibrium conditions dE T /dr = 0, 
dE T /dh = in normalized form as 

l( i +f)+l-i( i+ ¥) + ^-^ = ' (22a) 

| (i + x) + ~ iXJw + * x3 ~l 1{x) + F(a;) = °- (22b) 

To show the biasing effect of the top wall we write eq. (22a) as 

asffl + .[a + iJa]_ TO -o (23) 

where l(h/2) and Z(/i) are the material lengths as measured at the middle 
and at the top of the domain, respectively, and are denned by the 
obvious generalization of eq. (1). Equation (23) should be compared 
with the force equation for the standard cylindrical domain [Ref. 2, 
eq. (69)] 

\ + xff - F(x) = 0. (24) 

n 

Equations (22a) and (22b) may also be used to obtain expressions for 
the values of y a and ff required to obtain an equilibrium state charac- 
terized by specified values of l/h and x: 

t---is( 1 -|) + 5' w -" w (25a) 



1940 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1972 

a = -I I d + l + 1) - 2 (; + ?>« + (> + j)w. (25b) 

The stability or instability of a given equilibrium point is determined 
by examining the second derivatives of the energy. Calculating these 
derivatives and removing a factor of /3 -1 = kirh • ixM] we find 

* dr 2 a - 2 h + h ' + H ~ dx ' (26a) 



^^ = 2 A + 2M 1 + 2J + 2 XH - FW + 2"dT 



(26b) 



_ d 2 tf r 1 3 r/ . . „, v x 2 dF(afl /0 . . 

'if = i^» - 2 m + * F(x) - r & ' (26c) 

+ |tf-Jl W + S fW .lffl*. (2 6 d) 

Operations with these second derivatives are simplified if we define 

0(x) = -I(x) + sF(s) - I ^. (27) 

Using the relation 

- I(x) - F(x) = l(S.(x) + 3S 2 (z)) (28) 

we may write (?(x) in terms of the stability functions S„(x) and <S 2 (x) 
of Ref. 3: 

G(x) = §x(S (x) - S 2 (x)). (29) 

After consulting plots 3 of S (x) and S 2 (x), it is easily seen that G(x) ^ 
for x ^ 0. 

If we use eqs. (25a, b) to eliminate y w , H from eqs. (26a-d), we 
obtain expressions for the second derivatives of E T at equilibrium in 
terms of l/h and x only 

»(&1 - 41 & +*>->>- < 3 °>» 
<® -414-1) + §<**>■ (30c) 
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The derivatives are now in a form suitable for comparison with the 
stability conditions, eqs. (5-7). Inspection of eq. (30d) shows that 
P(d 2 E T /drl) eq , > for all positive values of l/h so that eq. (7) is satisfied 
for this system. Equation (6) is satisfied if 

i <?<*■>£ + sT- (31) 

Finally, we require 

JfafEA (d a E T \ _ ( d>E T Y 1 

which may be factored 7 into l/h > and 

i > g + 3)G W (l - , - ff. (33) 

However, it is easy to show that eqs. (31), (33) and the condition l/h > 
cannot be simultaneously satisfied. Thus, if y u = 0> the half-bubble 
model has no stability against collapse and/or run-through. 

While its result is negative, this stability calculation illustrates the 
method used in this paper for separating the l/h- and d//i-dcpcndence 
of the half-bubble model and for finding its stability region. We first 
use the equilibrium conditions to remove y w and H from the expressions 
for the second derivatives of E T . The resultant expressions are essen- 
tially polynomials in l/h with coefficients depending only upon x = d/h 
[cf. the right-hand sides of eqs. (30a-d) and (32)]. For fixed x, the roots 
of these polynomials provide us with the limits of the stability region 
(if one exists). In the next section, we generalize this method to include 
the effects of non-zero y M ■ 

TV. HALF-BUBBLE STABILITY WITH y u ^ 

We now generalize the equations of Section III to calculate the 
stability of the half-bubble model in the presence of a gradient in the 
saturation moment M, . The wall energy E w is unchanged from the 
expression given in eq. (10) and the interaction energy AE H only gains 
a new factor to become 



X < v - ; ~ [-5(2 - h) + |u(a)u(ft - 2)]}- (36) 
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AE„ = M.(l + ^)fl* jf " fft?) dp. (34) 

The expressions for the magnetostatic energy, however, become sub- 
stantially more complex than those given in eqs. (15) through (19). 
The magnetization distribution of eq. (9) is now replaced by 

M = Uf.(0)(l + ^)[1 - 2u(r b (<p) - r)]u{z)u(h - z). (35) 

The corresponding magnetostatic charge density is 
-V-M = -M.(0)[1 - 2u(r 6 (?) - r)] 

y_Mh 
I 

The first two 5-functions in this expression are the familiar platelet 
surface charge terms which occur in the stability theory of the standard 
cylindrical domain [Ref. 3, eq. (19)]. The third 5-function gives the 
correction to the charge on the upper surface, while the product of 
step-functions describes the uniform volume charge resulting from the 
gradient y M • 

The magnetostatic energy integral of eq. (14) has a quadratic de- 
pendence on the charge density — V • M. As a result, we have terms in 
AE M up to second order in y M . For convenience, we write 

AE M = tfff + yJB$ + faff (37) 

where E ( m is identical to the magnetostatic energy of the preceding 
section. The quantity y^Eff has two parts. The first arises from the 
interaction of the y M = charge distribution with the uniform volume 
charge and is zero by symmetry. The second, representing the inter- 
action of the y M = charge distribution with the excess charge on the 
upper surface, is easily evaluated. We find 

m TJT(l) yM"- 7J1(0) /OD\ 

yu&u — ~j &m • \po) 

The final term y M E l M may be written as 

faff = \ (u&Jm'M f~ r dr jT d<p f" r' dr' fj d<p'(kk' - 1) 

X [p -1 - 2V 1 £ dz[p 2 + (« - h) 2 ]* 

+ h~ 2 f dz f dz'[p 2 + (z - z') 2 A- (39) 

The p -1 term in the curly brackets corresponds to the self-energy of the 
excess charge on the upper platelet surface; the z-integral following it 
arises from the interaction of this excess charge with the uniform volume 
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charge, and the double-2-integral corresponds to the self-energy of the 
volume charge. After evaluating the z-integrals, we obtain the following 
integral expression for the second-order energy: 

Tid? = (^T^MliO) f rdr f d<p f r' dr' [" d<p'(kk' - 1) 

\ l I Jo Jo Jo Jo 

x (v - '-<£+*n . (40) 



Setting r 2 = 0, we find 
15 \ Z 



7« = T5 (^W^O))^ 3 ) 



f- {x + x*)F(x) - (3-x*)I(x)( Ml) 

The derivatives of the total energy .EV are now used to find the range 
of stability of the half-bubble model. Experimentally, this stable range 
would be mapped out by selecting materials characterized by given 
values of y w and y M and then varying 3 to find the stable region in the 
variables d/l and h/l. The analytic form of our equations, however, 
precludes following such a physical approach in the present calculation. 
Instead, we use an extension of the method followed in Section III. 
We use the equilibrium conditions eqs. (4a, b) to eliminate y w and ft 
from the second derivatives and then treat the second derivatives as 
polynomials in y = h/l, the coefficients of the polynomials being func- 
tions of x = d/h and y M . This manifestly unphysical approach gives 
the desired stable regions with a minimum of algebraic and numerical 
complexity. 

The equilibrium conditions corresponding to eqs. (4a, b) are con- 
veniently written in the following normalized forms: 

I ( a + 1) + * (1 + x) + ,4 + ad) _ ( x + f )f W 

+ I ^t)'\j ~ (1 + «VW + xl(x)j = 0, (42a) 

+ 3f [*» - I /(*)] - (^ff i m = 0. (42b) 

These equations may be solved to find the appropriate values of y w 
and ft for an equilibrium state characterized by given values of y u , x> 
and y = h/l. For brevity we write F(x) and I(x) as F and I respectively, 
and define a denominator D by 
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D = y[3 + 7„(l-f),]- (43) 



We then find 



y. = 0-'|(x - 2) + (y„x - 6F + 12x-I)y + r„(-10F + 22 X -'I)y* 
+ ly>„[£-(l + X >)F + (l + x )l]y<}, (44a) 

- 7 2 .[| (* + l)(;jr - (1 + «V + xl) + I (1 + z)/]*/ 3 }- (44b) 

The denominator D in both of these expressions is always greater than 
zero and finite in the region of interest, and thus contributes no sin- 
gularities or zeros to the calculation. 

The second derivatives of E T are developed in much the same fashion 
as in Section III. They are cast in dimensionless form by extracting a 
factor of /3 _1 = 4ir/i-47ril/^(0), and y w and H are then eliminated using 
eqs. (44a, b). The resultant expressions assume the form of polynomials 
in y after being multiplied by D. Because of the length of the expressions, 
we give only the final polynomial forms: 

--[( 4 -f)f-(HW3'> 

+4ff-( 2 -i+* ! )f 

+d-H>-&-4)'> 

+ (H + *>-(? +i » < 45a) 
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m 



d 2 E, 
dr„ dh 



= -|tt + 3)-[JG + f (l + ' 



y 



^,.i: 1 !.-i)s-(i-.)'+e-!)'> 






2/ dx- 



+ **' "S + 2 I 1 "" 2) dx 



_(|_|^ +[ . 



+ f|- 



DP 



dh 2 



= -111-:, 



i>/3 



d 2 E, 
drl 



i°- 2 f( l -i+f-;j» 

, . <iF , /15 a; ! \„ 23 .1 2 



[3 2 

a; a; 
6^~ 4 

+£[-s+(«— '+!>- (•-«•+ &> 

(45c) 

+ -[-(-!)f + te- i > + fe + i» 



+ 7jlf |(16 + 7z) — --^l-- + x - 



■[■ 



I^a-iZ_-_ 



+ f^ + i^-i- X V F + 1 (^ + 1 + f 



2/ dx 

(45d) 
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The problem of finding the ranges of parameters yielding stability 
of the half-bubble model now becomes one of numerical analysis. For 
given values of x and 7 W we find the zeros of the three polynomials 

r*M?E T \ (d 2 E T \ ( d 2 E T Y 1 

P L\ drl ) eq \ dh 2 U \dr dh) eq .]' 

MSI, - Mtl- 

A stable region corresponds to a range of y beginning and ending with 
one of these zeros and having the property that all three of the poly- 
nomials are positive for values of y on the interior of the interval. 
The results of this procedure are shown for y M = —0.025, —0.05, 
and —0.1 in Figs. 3-5. The stable regions are plotted in the x,l/y plane; 
1/z/ = Wh was used as a coordinate in order to facilitate comparison 
with the results of Thiele for the standard cylindrical domain (cf. Ref. 3, 
Fig. 3). In the interior of the stable region we have provided two families 
of curves which describe the behavior of a stable half-bubble. The 
families y w = constant and ff = constant are obtained from the poly- 
nomials of eqs. (44a, b) after the boundaries of the stability region have 
been found. Since h is not fixed in the half-bubble, the domains do not 
follow lines of l/h = constant as H is varied. Rather, they follow the 
curves y w = constant which characterize the given material. 

V. DISCUSSION 

We have explored the static stability of the half-bubble magnetic 
domain using a variational model. In spite of the simplicity of our 
model, the generality of some of the principles underlying our calcula- 
tion lead us to expect that our results have at least semi-quantitative 
validity (cf. Appendices A and B). 

Typical results of our stability calculation are shown in Figs. 3-5. 
The boundaries of each of the half-bubble stability regions shown in the 
figures roughly resemble two connected parabolic line segments, both 
parabolas opening to the right. The lower segment is determined by the 
zeros of the polynomial Df3(d 2 E T /drl) eq , and corresponds to the run-out 
instability mode. The upper segment arises from the zeros of 

This part of the boundary gives the limits of stability against run- 
through and collapse. Thus, the lines y w = constant on the right side 
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of the figure go from run-out at the low-/? end of the line to collapse 
on the high-i? end. Values of y w which are algebraically small enough 
to correspond to lines on the left end of the stable region follow a different 
pattern, going from collapse to run-through as H is decreased. Our 
stability charts then indicate two fundamentally different types of 
behavior as S is decreased: the half -bubbles may undergo run-through 
to assume the shape of a standard cylindrical domain (and then pre- 
sumably run out as H is decreased further), 8 or the half-bubble may run 
out first (presumably to undergo run-through as H is further lowered 
and give the usual demagnetized stripe pattern). These two types of 
behavior arc separated by a critical wall energy gradient y wc which is 
a function of y M , being approximately —0.027 and —0.04 for y M = 
—0.025 and —0.050, respectively. (For y M = —0.1, the problem is 
more complicated and will be discussed below.) 

Comparison of the half-bubble stability regions with those of the 
standard cylindrical domain (Ref. 3, Fig. 3) reveals one marked dif- 
ference between the two systems. In the half-bubble case, there is a 
minimum stable value of d/h for given y M . For y M = —0.025, —0.050, 
and —0.100, these minimum values are d/h = 0.39, 0.64, and 1.14, 
respectively. Long, narrow bubbles are then not attainable in the half- 
bubble system. 

It should be noted that all three of the stability charts correspond 
to 7 m < 0. While a general proof of the instability of half-bubbles with 
positive y M has not been formulated, numerical investigations show no 
regions of stability for y M > 0. As was shown to be the case for y M = 0, 
eqs. (5) and (6) could not be satisfied simultaneously f or y M > 0. This 
property of gradient-stabilized half-bubbles is actually an advantage 
for device work as it implies a useful type of stability. A half-bubble 
originating on the lower surface of a platelet, after being allowed to 
undergo run-through, cannot be squeezed down to a half-bubble residing 
on the upper surface as it is increased again. The domain is thus forced 
to return to its original state, the half-bubble state on the opposite 
surface being unstable. 

In addition to the stability limits imposed by the derivatives of the 
energy function E T , we place a lower limit on l/h by requiring a w {h) > 
and M,(h) > 0. While the motivation for the restriction <r w (h) > is 
physically obvious, the reason for the requirement M,(h) > stems 
more from the limitations of our model than from physical considerations. 
Magnetic materials may be compensated so that M,(z) goes linearly 
with z from positive to negative values. In the presence of a bias field, 
however, the magnetization in such a material would be aligned so that 



HALF-BUBBLES 1951 

its variation could no longer be described in terms of a linear variation 
with z. Such a compensated material would then not be susceptible to 
the present analysis and would presumably give rise to domains with 
different characteristics than those considered here. Taken together, 
the conditions c w {h) > and M,(h) > yield the restriction 

l/h > max (— y w , -y M ). (46) 

In practice, it may be necessary to require that l/h be even greater than 
the lower limit implied by eq. (46) since c w {z) > must be satisfied not 
only at the top of the half-bubble (z = h), but throughout the region 
of the platelet above the half-bubble (z > h). (By a simple extension 
of Appendix A, we need not require M t (z) > in the region above the 
half-bubble.) These restrictions indicate that our calculation is not valid 
for part of the lower end of the stability region obtained for y M = 
—0.100 (cf. Fig. 5). As a consequence, the critical wall energy gradient 
7„ c discussed above cannot be defined in this case. 

The normalized bias fields H required to stabilize half-bubbles are 
somewhat smaller than those needed for standard cylindrical domains 
as a result of the self-biasing effect provided by the upper closure wall 
(cf. Section III). Without going to unreasonably large values of y w , 
cases are found which are stable with H = or even H < (H parallel 
to the magnetization inside the domain). In addition, bias margins for 
the half-bubble are of the order of those for the standard cylindrical 
domain, typically falling in the range 0.04 < A/? < 0.08. Thus if 
4tM,(o) is 200 gauss, bias field margins will be of the order 8 to 16 gauss. 

To demonstrate the use of the stability charts in obtaining numerical 
information, we consider the case y M = —0.100 and y w = +0.025, for 
which a half-bubble is stable in zero bias field. The variations in a w 
and M, between the bottom and the top of the half-bubble are easily 
computed after obtaining a value of l/h from the stability chart. In the 
present case, l/h = 0.152 corresponds to zero bias. Evaluating y w h/l 
and yyfh/l, we find that a w increases by 16.4 percent and M, decreases 
by 65.8 percent between the bottom and the top of the half-bubble. 

The collapse and run-out fields for y M = —0.100 and y w = 0.025 are 
found by interpolating between the H = constant curves in Fig. 5. 
We find jff.,,11. = 0.017 and H Taa . ^ -0.023 giving a bias margin of 
AH = 0.040. The size variations of the half-bubble over its stable range 
are best measured in terms of the z = material length defined by 
eq. (1). For any operating point in the stable range, h/l = (l/h) -1 can 
be found from the vertical coordinate of the point while d/l = (d/fyil/h)' 1 
can be found as the ratio of the horizontal and vertical coordinates. For 



1952 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1972 

7M = -0.100 and 7* = +0.025 we find (h/l) ooU . = 5.38, (d/l) eoli . = 
10.85, (h/l) Tun . = 7.27, and (d/l) rm . = 30.91. Thus, while d increases 
by a factor of three between collapse and run-out, h increases by only 
35 percent over this same range. Since the curves S = constant and 
y w = constant are nearly parallel at the collapse end of the stable range, 
much of this change in h occurs with a very -small change in S near 
collapse. Thus, the relative variation in h may be decreased considerably 
by sacrificing a small fraction of the bias margin. The relative smallness 
of this variation in h is particularly important if half-bubbles on the 
lower platelet surface arc to be propagated by permalloy circuits at the 
upper surface, as it is to be expected that there is an optimum circuit- 
domain separation for half-bubbles much as there is for standard 
cylindrical domains. 9 

In conclusion, we have shown that the half-bubble has a number of 
properties which are attractive for device applications. Its stability in 
the absence of external biasing represents a considerable advantage over 
the externally biased standard cylindrical domain. While some com- 
plexity is added by the requirement of gradient materials, the materials 
problem is simplified on another front by the fact that material thickness 
is unimportant for half-bubbles. Indeed, the half-bubble stabilizes on 
only one surface; if propagated from the same surface, irregularities in 
the second surface are of no consequence. These attractive properties 
of the half-bubble and its possible utility in avoiding the undesirable 
properties of hard-bubbles may more than compensate for the added 
difficulty in growing half-bubble materials. 
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APPENDIX A 

Magnetostatic Energy Invariance 

Demagnetization fields equivalent to those of the standard cylindrical 
domain may be conceptually created through the use of two neutral 
subsystems of charges. The first consists of two parallel planes with 
magnetostatic charge density ±M, and the second consists of two 
disks with charge density T2M", placed on a common axis just inside 
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Fig. 6 — Source-charge distribution for the demagnetization field of: (a) The stand- 
ard cylindrical domain, (b) The half-bubble model. 



the first subsystem (Fig. 6a). The distances between the upper plane 
and the upper disk, and between the lower plane and the lower disk 
are allowed to approach zero. The demagnetization field of the cylin- 
drical domain is the superposition of the fields of the two isolated 
subsystems: 

H d = H dl + H d2 . (47) 

The magnetostatic energy of the domain (including the energy of the 
uniformly magnetized platelet) is the sum of the self-energies of sub- 
systems 1 and 2 and their mutual interaction energy: 



E t 



Uv 



= hf» 



-kl*» dv + hl 



H 2 d2 dV + 



4tt J 



H< t -H M dV. 



(48a) 



(48b) 



The transition from the magnetostatic field configuration of the 
standard cylindrical domain to that of the half-bubble model is achieved 
by moving the upper plane of charge upward (Fig. 6b). The field strength 
of subsystem one is unaltered by this change, while the self-energy of the 
subsystem increases linearly with the distance moved. The self-energy 
of subsystem two is, of course, unchanged, but remarkably the mutual 
interaction energy of the two subsystems is unchanged also. The latter 
invariance follows from the invariance of the field strength of subsystem 
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one and of the charge-neutrality of subsystem two (so that all field 
lines originating in the isolated subsystem two also terminate there). 
From these two properties it is easily seen that any lines of R d2 extending 
above the upper disk give a zero contribution to the mutual interaction 
energy of the two subsystems. If we now observe that the magnetostatic 
energy of introducing a domain into a uniformly magnetized platelet 
is just the sum of the self-energy of subsystem two and the mutual 
interaction energy, we see that this differential magnetostatic energy is 
independent of the position of the upper surface of the platelet (provided 
only that the domain remains inside of the platelet). 

It is easily seen that this proof may be generalized to the cases in 
which M, is a function of z or in which the domain assumes a more 
general form. 

APPENDIX B 

Energy of Half-Bubble Under Elliptical Distortions 

The condition (dE T /dr 2 ) Tl=0 = in the absence of in-plane anisotropy 
is satisfied not only by our simple half-bubble model but by any convex 
domain having radial symmetry at equilibrium. For example, consider 
the generalized half-bubble shape 

r b (<p, z) = r (z) + r 2 (z) cos 2«p (49) 

where r 2 (z) « r„(z) and r (h) = (i.e., the half-bubble is closed at the 
top). 

If there is no in-plane anisotropy, then 

E T (r 2 ) = E T (-r 2 ) (50) 

since the shape function r b (<p, z) = r„(z) — r 2 (z) cos 2<p transforms into 
that of eq. (49) under rotation. From eq. (50) it follows that 

dE ^ ^0. (51) 

Consequently, cross-derivatives such as (d 2 E T /dr 2 dh) r ,. are also 
zero for the general domain shape given in eq. (49) . 
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